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A non-perturbative quantization of a paraxial electromagnetic field is achieved via a gener- 
alized dispersion relation imposed on the longitudinal and the transverse components of the 
photon wave vector. This new theoretical formalism yields a seamless transition between the 
paraxial- and the Maxwell-equation solutions. This obviates the need to introduce either ad hoc or 
perturbatively-defined field operators. Moreover, our (exact) formalism remains valid beyond the 
quasi- monochromatic paraxial limit. 
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PACS numbers: 03.65.Db, 03.70.+k, 41.85.-p 

In this Letter we seek an answer to the question: What 
is the quantum-mechanical state of a photon in a beam of 
light whose propagation can be classically described by a 
paraxial wave equation? The problem of quantum prop- 
agation of paraxial fields was considered by several au- 
thors in the past Q] . In these contributions the approach 
was either ad hoc or approximated. A notable excep- 
tion was the work of Deutsch and Garrison where the 
authors developed a perturbative quantization scheme 
0. However, although their theory formally solves, or- 
der by order in a perturbation expansion, the problem of 
paraxial quantum propagation, it suffers from two main 
limitations. First, it does not provide any clear and eas- 
ily manageable formula for the paraxial quantum modes. 
Second, it requires the quasi-monochromatic approxima- 
tion which is unsuitable for, e.g., the description of the 
very quantum phenomenon of propagation of broad-band 
entangled photons provided by spontaneous parametric 
down conversion (SPDC). The increasing importance, for 
quantum information in general 0] and quantum cryp- 
tography £| in particular, of a proper description of 
paraxial propagation of SPDC entangled photons, calls 
for an exact (namely, non-perturbative) theory of parax- 
ial quantum fields. 

The aim of this Letter is to introduce a non- 
perturbative method to quantize a paraxial electromag- 
netic field. Our scheme is conceptually simple: we begin 
by considering the quantized transverse electromagnetic 
vector potential calculated at the (arbitrary) initial time 
t = 0. Then we select from the wave-vector space only 
those field configurations which are exact solutions of the 
paraxial wave equation, and let them evolve in time ac- 
cording to the d Alembert wave equation. This procedure 
automatically ensures the validity of both the canonical 
commutation relations and the transversality conditions 
for the fields. We stress that our theory accounts for 
the axial propagation of vector fields with any spectral 
and spatial bandwidth and reproduces the well known 
paraxial results in the limit of narrow beam-like fields. 

To begin with, let us consider the vector potential op- 
erator A(r,t) = A' + '(r,t) + A(~)(r,i) in the Coulomb 
gauge which can be written in the plane-wave basis as 



A( + '(r,t) = / d 3 k 
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A=l 
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(A) (k)a A (k)exp (ik • r ■ 



ic|k|f), 



(1) 

and A( ) (r, t) is the Hermitian conjugate of A^ + ) (r , t) . It 
clearly satisfies the d'Alembert wave equation DA(r, t) = 
0. The two unit polarization vectors e( A )(k),(A = 1,2) 
are transverse e^(k) • k = and mutually orthogo- 
nal eW(k) • e( 2 )(k) = 0. Moreover, the annihilation 
and creation operators satisfy the canonical commuta- 
tion rules [a A (k),a A ,(k')] = <S A A'<5 (3) (k - k'). Since the 
fields considered in this paper are mainly beams which 
propagate close to the z-direction, we find convenient 
to introduce a finite quantization length L along the z 
axis with discrete wave vector Cartesian z-components 
k z — > ( n = jln, (n = 0,1,...). The choice n > 
implies that only the parts of the field which propa- 
gate in the positive z-direction are included in Eq. JQ. 
The integral with respect to dk z in Eq. <QJ is then re- 



placed by a sum |fj, [7| / dk x dk y dk z d 2 k T , 
where ky = (k x , k y ) and d 2 k^ = dk x dk y . Moreover, we 
scale the annihilation and creation operators by defining 
a\(k) = yj L/ (2ir)a\{kT , n), in such a way that 

{a x (k T ,n),al(k' T ,n')} = 6 xx ,5 nn ,6W(k T - -k' T ), (2) 



where the Kroncckcr symbol replaces the delta function 
according to 5{k z — k' z ) — > ^S nn >. 

A paraxial field is usually expressed as an envelope 
field modulating a carrier plane wave with wave-vector 
k and angular frequency ojq = c\k \. Without lack of 
generality we assume k = k z and we choose k > so 
that the carrier plane wave propagates in the positive z 
direction. By using the trivial identity 

C f^ /C A , exp(ifc z - kj t) 

1 = 7T / dfc 771 : — 77, (3) 

Z7T J Q exp(ifcoz — iu>ot) 

where C is an arbitrary length, we can rewrite the vector 
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potential operator as 
P(r,i) = A(+»(r,t) x 1 

dko exp(ifco^ — iuJot)^f(r, t), 



2tt 



where we have introduced the envelope field 



*(r,t) = E/ ^ 



1/2 



(4) 



(5) 
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A=l,2 

x exp [i(k — fc z) • r — i(c|k| — u>o)t], 



where k = ky + z£„. Now, the key idea is to find a 
subspace of the three-dimensional wave vector k-space 
where the initial time envelope field *(r) = *(r,< = 0) 
satisfies the paraxial wave equation Q 

d 2 4>(r) | d 2 *(r) | g.^ A* (r) _ Q (g) 
9x 2 ch/ 2 <9z 

If we substitute from Eq. (JSJ into Eq. © we obtain 

( n /k = l-\k T \ 2 /2k 2 . (7) 

This generalized dispersion relation plays a key role 
through this paper. It defines a two-dimensional domain 
in the k-space where both d'Alembert wave equation and 
paraxial wave equation are satisfied. For sake of clarity 
we write q = k^, q — | q| and define the dimensionless 
parameter = q/(^/2k^). If we denote with 9 6 [0,7r/2] 
the angle between the wave vector k and the axis z, 
then g/k • z = tan# and from Eq. JJJ), it follows that 
= -cot8 + yj{2 + cot 2 9). This relation is exact; 
however, we can gain some insight if we consider it in the 
limit 9 < 1 where $^2 ~ 9 - 9 3 /6 + 0(6> 5 ). This equa- 
tions shows that "3^/2 = q/ko is approximatively equal to 
the divergence angle of a Gaussian beam Moreover, 
by comparing $ with Eq. (2.8) by Deutsch and Garri- 
son 0, one recognizes z?\/2 < 1 as their perturbative 
expansion parameter. However, in our case the only con- 
straint is § < 1, as follows from Eq. J7J) and the condition 

Cn>0. 

From Eq. Q) it readily follows that the exponential 
factor in Eq. J5| can be written 



exp [i(k — fc z) • r — i(c|k| — uJo)t] 

= exp(iq • x - i$ 2 k z) exp[-iw *( \A + i? 4 - 1)], 



(8) 



where x = (x, y). Clearly Eq. (JTJ affects also the value of 
the polarization unit vectors in Eq. (JSJ). To see this, we 
first write the total wave vector k in terms of q and ko as 
k = qg+zfco(l— # 2 ), where q = q/|q|, then we arbitrarily 
choose (we always have this freedom) (q, d) = z x q. 
The remaining unit vector is then uniquely fixed by 
the cyclic relation cx x k to the value 

e«(q,tf) = [q(l-^ 2 )-z ? 9V2]/(l + ^ 4 ) 1 / 2 . (9) 
It is now possible to write explicitly the envelope field 
^(r) restricted to the k-subspace defined by the disper- 
sion relation Eq. (JTJ). From the definition of and 
Eq. Q, it follows that we must select from the sum 
over n in Eq. (JSJ) only those terms corresponding to 
n = n($) ee (l - tf 2 )] Ip , where "IP" stands for In- 
teger Part. This objective can be achieved by replacing 



in Eq. © E, 



where <L 



= 1 for 



n = n("&) and 8 n n ^ — otherwise. Finally, we can 
write from Eq. JSJ 



*(r) 



h 



1/2 



r 2 £ w (l + ^ 4 ) 1 /2 J L / 
£ e^\q^)a x (q,n)5 nM ^ ( 10 ) 



A=l,2 



x exp(iq ■ x — i-d koz), 

where e' A ) (q, $) are given by Eq. Q and previous for- 
mulae. We have substituted everyplace in Eq. lfTT)|l by 
Eq. 0; this operation is permitted by the presence of 
^n,n(i?) within the sum. The only exception to these sub- 
stitutions is represented by the operators a\ (q, n) which, 
for the moment, are left unchanged. 

At this point, we note that since the restriction to the 
paraxial k-subspace has already been achieved via Eq. 
(|10ll . it is possible to make a step backward from the 
discrete momentum C, n to the continuous frequency lj: 
c£ n = 2 -j^n — > uj. The other required replacements are 

£„ ~> 2^/ dw ' «A(q,n) = y^pa A (q,w) and 
2irr 

S n ,n W = —5[lu- ck (1 - $ 2 )] , (11) 
in such a way that 

[a A (q, u), &{, (q', u')] = 8».>S(U> - lo')5^ (q - q'). (12) 
Equation Q can be then written as 



A' + '(r./l = (§ 



2-n/C 



dfco exp(ifco2 — it^oi) / dw / d q 



o 

x e(A) (q> ^(q, oj)5 [w - cfc (l - $ 2 )] exp(iq • x - id 2 k Q z) exp[-iw i(\/l + ?? 4 - 1)] 

A=l,2 



n 



167r3e fc (l + 1 ? 4 ) 1 / 2 



1/2 



(13) 
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where the first term C/L is a dimensionless constant fac- 
tor which can be eliminated by renormalizing A( + ) (r, i) 
in the end of the calculations; therefore we leave it out 
from our formulas. This is not yet our final expression 
since we can perform explicitly the integration with re- 
spect to kg by using the well known formula for a delta 
of a function which gives us 

5 [u - ck (i - tf 2 )] - s {ck a - n )/(i + e 2 ), (w) 



where 

6 = v /2g 2 c 2 /[o; + \A> 2 + 2q 2 c 2 }, (15) 

flo/c = q/(®\2) and only the term corresponding to 
fc > has been retained. Finally, we can write 



A i+) (r,t) = J dcjexp[-iw(i- z/c 



167T 3 £oCW 



1/2 



J d2c t Yl ^ (A) (q^^^)aA(q,w)exp ^iq-x-ij^^, 



r 



where we have defined the slowly varying polarization 
vectors 

1/4 

£( A >(q, W ,M)^e( A )(q,e) 



(i + e 4 )(i + e 2 ) z 



xexp 



(lo - fWi + e 4 ) t-i(u)-n ) z/c 



(17) 

and the unit vectors e( A )(q, 0) are given by Eq. © and 
previous formulae with •& —> 0. 

Equation l|16l) is the first main result of this work. It 
represents a vector potential field operator which is a 
bona fide transverse field obeying the d'Alembert wave 
equation for any time t > 0, whose corresponding enve- 
lope field satisfies the paraxial wave equation |J§} at t = 0. 
We stress that this expression is exact, no approximations 
were made. 

Now, we are ready to address the problem of build- 
ing the quantum mechanical state describing a photon 
in a paraxial beam. To this end, we first note that un- 
til now the creation operators fit(q, u>) have passed un- 
touched through all our operations. In fact, they still 
satisfy the canonical commutation relations Eq. (|12J) . 
However, their form is not the most suitable one to deal 
with paraxial fields; therefore we introduce the Fourier- 
transformed operators 0, fiof 



a A (x,w) 



— J d 2 qa A (q,cj)exp(iq- x), 



(18) 



such that 



[a A (x,c),ai,(x',c')] =<5aa^(^-c')^ 2) (x-x'). (19) 
ifT^I in Eq. JUJ we obtain, after 



If we substitute Eq. 
some algebra, 



A (+) (r,t) = j / dau 



A=l 



-iu)(t— z/c) 



(Aireocw/h) 1 /' 



i(A). 

A {x,z,uj,t), 



(20) 



where r = (x, z) and we have introduced the exact slowly 
varying photon annihilation vector operators 

A (A) (x,z,w,i) = J d 2 x'^ A )( x , z,x' ,w,t)ax(x! ,w). 

(21) 

Moreover, in Eq. I|21|l we have defined the Maxwell- 
paraxial (MP) slowly varying modes 



^ (A) (x,z,x>,t) 



^5 J d 2 q^(q,^^i) 



(27T) 

xexp 



Iq ' (x - X) -W 



(22) 

Equation i|22[l displays the second main result of this pa- 
per. It describes the field in the plane z at time t due 
to a point source with frequency u> located at x' in the 
transverse plane z = 0. As its shape clearly suggests, 
^ (A) (x,z,x',w, t) is the quantum analog of the classical 
Huygens-Fresnel diffracted field. This may be seen more 
clearly by writing Eq. i|22|) in the narrow-beam limit 
which is achieved by restricting the transverse momen- 
tum integral to the domain C w = {q : q <C u)/c}. It 
is easy to see that within this domain ~ <C 1, 

n ~ uj and £ (A) (q,w, z,t) ~ e< A )(q), where eW(q) = 
q + O(0) and e( 2 )(q) = z x q, are the zeroth order polar- 
ization unit vectors. In this limit Eq. 12211 reduces to the 
well known paraxial Green's function 7 > ^(x, z, x', u>) || 



7> (A) (x,z,x» 



d 2 q 



,(A) 



(q)e 



i q-(x-x') 



(23) 

here generalized to vector fields. Now, it is straightfor- 
ward to show that at a fixed time t, in each transverse 
plane z, the MP functions J-^ 1 (x, z, x', u>, t) are quasi- 
orthogonal 



4 



J d^^\^ z y,u,t)-^\^zy,^t) = s x ,J ( ( i + e^) (fl ' ex pfa-( x '- x ")]< ( 24 ) 



that is, the right side of this equation approaches 
Sx^S^ (x' — x") in the narrow-beam limit u/(qc) — > oo, 
as expected from the orthogonality of the classical parax- 
ial Green's functions for free space propagation. This 
result was already found by Visser and Nienhuis 0] 
who suggested to interpret the Fourier-transformed cre- 
ation operator a\(x.,u>) as the operator which creates 
at t = a photon with polarization A in the paraxial 
mode exp(kjz/c)'P ( - A ' ) (x, z, x', u). More generally, the in- 
terpretation of ^^(x, z, x', cj, t) as single-photon wave 
function can be put on a rigorous basis by introducing 
the "transverse-position" states |x, cj,A) = a^(x, w)|0). 
Then it readily follows that (0|A^ + '(x', z, i)|x, lu, A) oc 
F (X) (x', z, x, u, t) exp[-iw(t - z/c)} . 

Now, the definition of |x, w, A) makes possible to as- 
sociate to any single-photon state 1^), its corresponding 
Maxwell-paraxial wave function tp\(x., u>) = (x, w, X\ip). 
Then, for example, the MP wave function associated 
to the plane wave state \q,ui, A) = a^(q,o;)|0), is sim- 
ply given by the Fourier relation (x, a/, A'|q, lj, A) = 
cxp(iq-x) g^fi^ _ ^/y More generally, for a given com- 
plete set of orthogonal transverse functions ip nm (x, uj) as, 
e.g., the Hermite- or the Laguerre-Gaussian beams j9(, it 
is possible to build the corresponding MP single-photon 
state as 

\n,m,uj,X) = J d 2 x^, lm (x, w)|x,w, A), (25) 

where , X'\n,m,uj, X) = Vw(x> u)5\\<5(uj - w'). 

This equation is our third and final main result: it 
represents the exact quantum-mechanical state of each 
of the photons in the classical paraxial beam ■0„ m (x, lo). 

In conclusion, a non-perturbative quantization scheme 
for electromagnetic paraxial fields has been introduced. 
It relies on the fact that it is possible to select some ini- 
tial field configurations which are exact solutions of the 
paraxial wave equation. These configurations are then 
evolved at later times with the d'Alembert wave equa- 
tion. In this way we were able to find explicit and man- 



I 

ageable expressions for the exact field in both momen- 
tum [Eq. l(TH|l ] and in position [Eqs. (|20I22|) ] represen- 
tation. Moreover we gave an unambiguous definition for 
Maxwell-paraxial quantum states [Eq. (|25|l ]. This quan- 
tization method also suggests that beyond the paraxial 
case a whole class of other exact solutions of both parax- 
ial and d'Alembert wave equations could be found. 
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